Gluon distribution in proton at soft and hard pp collisions 
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^ I Abstract 

O I We analyze the inclusive spectra of hadrons produced in pp collisions at high 

energies in the mid-rapidity region within the soft QCD and perturbative QCD 
assuming the possible creation of the soft gluons at low intrinsic transverse momenta 

C/^ . kt- From the best description of the LHC data we found the parametrization of the 

! unintegrated gluon distribution which at low kt is different from the one obtained 

^ within the perturbative QCD. 

Qh' 1. Introduction 

^ . 

I— I. As is well known, hard processes involving incoming protons, such as deep-inelastic 

^ \ lepton-proton scattering (DIS), are described using the scale-dependent PDFs. A 

>• I distribution like this is usually calculated as a function of the Bjorken variable 

^ ' X and the square of the four- momentum transfer = —Q^, integrated over the 

Tij- ■ parton transverse momentum kt. However, for semi-inclusive processes, such as 

^ 1 . inclusive jet production in DIS, electroweak boson production [1], etc., the parton 

^ I distributions unintegrated over the transverse momentum kt are more appropriate. 

y—i [ The theoretical analysis of the unintegrated quark and gluon PDFs can be found, 

^ ' for example, in [2]-[5]. According to [3], the gluon distribution function g{kt) at 

^ ■ fixed has the very interesting behaviour at small x < 0.01, it increases very 

^ . fast starting from almost zero values of kt. In contrast, the quark distribution 

q{kt) is almost constant in the whole region of kt up to kt ~ 100 GeV/c and much 
smaller than g{x). This parametrization of the PDFs was obtained in [3] within 
the leading order (LO) and next-to-leading order of QCD (NLO) at large from 
the known (DGLAP-evolved [6]) parton densities determined from the global data 
analysis. At small values of the nonperturbative effects should be included to 
get the PDFs. The nonperturbative effects can arise from the complex structure 
of the QCD vacuum. For example, within the instanton approach [4] the very fast 
increase of the unintegrated gluon distribution function at < fct < 0.5 GeV/c 
and = 1 (GeV/c)^ is also shown. These results stimulated us to assume, that 
the unintegrated gluon distribution in the proton can be included by analyzing also 
the soft hadron production in pp collisions. In this paper we analyze the inclusive 
spectra of the hadrons produced in pp collisions at LHC energies in the mid-rapidity 
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region including the possible creation of soft gluons in the proton. We estimate 
the unintegrated gluon distribution function (UGDF) at low intrinsic transverse 
momenta kt < 1.5 — 1.6 GeV/c and its parameters extract from the best description 
of the LHC data at low transverse momenta pt of the produced hadrons. We also 
show that our UGDF is similar to the UGDF obtained in [2, 5] at large values of kt. 



2. Inclusive spectra of hadrons in pp collisions 
2..1 Quark-gluon string model (QGSM) 

Let us analyze the hadron production in pp collisions within the QGSM [7] or the 
dual parton model (DPM) [8] including the transverse motion of quarks and diquarks 
in colliding protons [9, 10]. As is known, the cylinder-type graphs presented in 
Fig.l make the main contribution to this process. The left diagram of Fig.l, the 
so-called one-cylinder graph, corresponds to the cut one-pomeron graph and the 
right diagram of Fig.l, the so-called multi-cylinder graph, corresponds to the cut 
n-pomeron graphs. The general form for the invariant inclusive hadron spectrum 




Figure 1: The one-cylinder graph (left) and the multi-cylinder graph (right) 
within the QGSM is the following [7]: 

p{x,pt) = E— = an{s)(f)n{x,Pt) , (1) 
P n=l 

where E, p are the energy and three-momentum of the produced hadron h in the l.s. 
of colliding protons respectively; x,pt are the Feynman variable and the transverse 
momentum of h; an is the cross section for production of the n-pomeron chain (or 
2n quark-antiquark strings) decaying into hadrons, calculated within the "eikonal 
approximation" [11] the function 0„(x,Pf) has the following form [9]: 

(j)n{x,pt)=f dxi f dx2ipn{x,pt;xi,X2) , (2) 

where 



1pri{x,Pt; Xi,X2) = 
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F^^\x+,Pt;x,)F^::\x_,p,;x,)/F(:\0,p,) + 
+F^:\x+,pt;x,)F^^\x^,Pt;x2)/F^^\0,Pt) + 
2{n - l)F(^\x+,pt;x,)F^:\x^,pt;x2)/Fi:\0,pt) . (3) 

and x± = 0.5(^x2 + xf ±x),Xt^ '^\/{ml + Pt)/s, 

I d%f!r\x±, kt)Gr^H [^/^t;Pt)^ , (4) 

F^''\0,Pt) = 

[' dx'd%f^^\x', kt)G,^H{0:Pt) = Gr^h{0:Pt) • (5) 

Here r means the flavour of tlie valence (or sea) quark or diquark, f^^\x'^ kt) is the 
quark distribution function depending on the longitudinal momentum fraction x' 
and the transverse momentum kf in the n-pomeron chain; Gr^h = zDT^h{z, kt]Pt), 
Dr^h{z,kt;pt) is the fragmentation function of a quark (antiquark) or diquark of 
flavour r into a hadron h. 



2.. 2 Hadron production in the mid-rapidity region 

. According to the Abramovskiy-Gribov-Kancheli cutting rules (AGK) [12], at mid- 
rapidity only Mueller-Kanchcli type diagrams contribute to the inclusive spectrum 
of hadrons. It was shown in [9, 13] that at a; = the function F!jP'\x± ~ 0,pt, 0:1^2) 
does not depend on n. Therefore, the function 0n(O,Pt) becomes proportional to n. 
In accordance with [13], the inclusive spectrum corresponding to the contribution 
of quarks at x = has the following simple form: 

00 

P<z(0,Pt) = 0g(O,pt) nan{s), 

n=l 

g(s/so)^U^,Pt), (6) 

where 0^(0, pj) is the so-called vertex function in the Mueller-Kancheli type diagram 
which is determined by the density distribution of the particles produced in the 
central region; A = 0.12, g=21 mb and and is the nondiffractive cross section due 
to the exchange of any number cut-pomerons. 

According to Eqs.(2-3), the sea quarks contribute to (;6„ and the inclusive spec- 
trum at n > 2. Assuming possible creation of soft gluons in the proton, which 
are split into qq pairs and should vanish at the zero intrinsic transverse momentum 
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{kt ~ 0), one can obtain the following energy dependence for the gluon contribution 
Pg{0,Pt) to the inclusive spectrum [13]: 

OD 

Pg{x = 0,pt) = J2 ^n{s)(j)U^,Pt) = 
n=2 

Finally, we can present the inclusive spectrum at x ~ in the following form: 

p{Pt) = Pq{x = 0,pt) + pg{x = 0,pt) . (8) 

The functions 0g(O,pt) and (j)g{0,pt) are fixed from the fit of the experimental data 
on the inclusive spectra of the charged hadrons produced in the central pp collisions 
at different initial energies. So, we found the following parametrization for 0^(0, pt) 
and (f)g{0,pt) [13]: 

0?(O,Pt) ^ Aqexp{-bqpt) 
^g{Q,Pt) = Ag^/p'texp(-bgPt), (9) 

charged particles at y = 0: where Ag = 4.78 ± 0.16 (GeV/c)-^ bq = 7.24 ± 0.11 
(Gev/c)-^ and Ag = 1.42 ± 0.05 (GeV/c)-^; bg = 3.46 ± 0.02 (GeV/c)-i. 

2.. 3 Gluon distribution in proton within GBW model 

The conventional QGSM or DPM model does not include the distribution of gluons 
in the proton. However, as is well known, at large transverse momenta pt of hadrons 
the gluons in the proton play very important role in description of the experimental 
data. Therefore, one can assume that the contribution of the gluon distribution in 
the proton to the inclusive spectrum of the produced hadrons slowly appears when 
Pt increases and it will be sizable at high values of pt. This assumption is also 
confirmed by the increase of the UGDF in the proton at a; ~ as a function of the 
internal transverse momentum kt when kt grows [2]- [5]. 

According to Refs.[2, 5], the UGDF, as a function of kt at some value of Qo and 
low X is presented in the following form: 

xg{x, kt, Qo) = CoRl{x)k^t exp , (10) 

where Cq = 3(To/ (4vr^Q;s(Qo))) Ro{x) =GcY^^ {x / xq)'^/'^ , xq, A and (Tq are defined in 
[2, 5]; as{Qo) is the QCD coupling constant. To get the UGDF at low kt we assume 
the possible creation of the soft gluons in the proton which appears at nonzero kt- 
We calculated the gluon contribution 4>g{0,pt) entering into Eq.(7) as the cut graph 
of the one-pomeron exchange in the gluon-gluon interaction (Fig.2, right) using the 
splitting of the gluons into the qq pair. Then the calculation was made in a way 
similar to the calculation of the sea quark contribution to the inclusive spectrum 
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Figure 2: The one-pomeron exchange graph between two gluons in the elastic pp 
scattering (left) and the corresponding cut graph (right). 



within the QGSM, see Eq.(3) at n = 2. Calculating the diagram of Fig.2 (right) we 
assumed the following form for the xg{x, kt, Qo)'- 

xg{x, kt, Qo) = CoCi{l - xf^ x 
{Rl{x)k^ + C2{Ro{x)ktr)exp{-Ro{x)kt) , (11) 

where C2, a, bg and A are the parameters which were found from the best description 
of data. As it will be shown in the next section, the parametrization given by Eq.(ll) 
results in the satisfactory description of the inclusive spectra of hadrons produced 
in pp collisions at the mid-rapidity region and low transverse momenta pt < 1.5 
GeV/c. At higher values of pt the gluon contribution (pg{0,pt) given by Eq.(9) is 
decreasing very fast and, in fact, does not contribute to the hadron spectrum at 
large pt- Therefore, one can find the parametrization for the UGDF that is very 
close to the GBW parametrization [2, 5] at large kf and coincides with the UGDF 
given by Eq.(ll) at low kt- It is the following: 

xg{x, kt, Qo) = CoCsil - x)*« x 
{Rl{x)k^ + C2{Ro{x)kt)'') exp {-Ro{x)kt - d{Ro{x)kt)^) , (12) 

The coefficients Ci,C3 are found from the normalization 

xg{x,Ql)= dk^xg{x,kt,Q^) (13) 
Jo 

3. Results and discussion 

The results of our calculations of the charged hadron inclusive spectrum 
{l/Nev)d^N/dydpt, where N^v is the total number of the inelastic events, in the 
mid-rapidity at \/i =7 TeV are presented in Fig. (3). The solid curve corresponds to 
the quark contribution pg{x = 0,pt) given by Eq.(6), the long-dash line is the gluon 
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contribution pg{x = 0,pt) (7) to the inclusive spectrum; the dotted curve corresponds 
to the sum of both contributions, see Eqs.(8). One can see that the conventional 
quark contribution pq{y = 0,pt) is able to describe the data up to pt <1 GeV/c, 
whereas the inclusion of the gluon contribution allows us to extend the range of good 
description up to 2 GeV/c. At larger values of pt the contribution of hard processes 
is not negligible. Therefore, we calculated the inclusive spectrum at pt > 2 GcV/c 
within the leading order perturbativc QCD (LO PQCD) [15]- [22], see the details in 
[13]. The short-dash curve in Fig. 3 (top) corresponds to our LO PQCD. Prom the 
best fit of the inclusive pt-spectra of the charged hadrons produced in pp collisions 
at X ~ presented in Fig.3 (bottom) by the short dash line we found the following 
parameter values entering into the form for xg{x, kt, Qo) given by Eq.(ll): 

a = 0.7; C2 ~ 2.3; A = 0.22; bg = 12; d = 0.2; C3 = 0.3295 

and 

Ci = ^(l + C2r(a + 2)/r(4))-^ . 

In Pig.4 the UGDF xg{x, kt, Qq)/Cq as a function of kt at x = xq and Qo = l.GeV/c 
is presented. The solid line (GBW UGDF) is the GBW parametrization [2, 5] given 
by Eq.(lO) divided by Cq; the dashed curve (UGDFl) corresponds to the UGDF 
given by Eq.(ll) and the dotted hne (UGDF2) corresponds to the UGDF given by 
Eq.(12) which is very close to the GBW UGDF at h > 1.5 GeV/c. It is seen from 
Fig. 4 that the parametrization UGDF2 is different from the GBW UGDF mainly 
at kt<lA GeV/c. 

4. Conclusion 

We assume that the contribution of the gluon distribution in the proton to the 
inclusive spectrum of the produced hadrons slowly appears when pt increases and it 
will be sizable at high values of p^. This assumption is also confirmed by the increase 
of the unintegratcd gluon distribution in the proton at a; ~ as a function of the 
internal transverse momentum kt when kt grows [3, 4, 2, 5]. 

Therefore, to illustrate this hypothesis we fit the experimental data on the in- 
clusive spectra of charged particles produced in the central pp collisions at energies 
larger than the ISR starting by the sum of the quark contribution pg given by Eq.(6) 
and the gluon contribution pg (see Eq.(7)). The parameters of this fit do not depend 
on the initial energy in that energy interval. From the best fit of the LHC data on the 
inclusive spectra of the charged hadrons produced in the mid-rapidity pp collisions 
at low Pt we found a new parametrization of the unintegrated gluon distribution 
in the proton at small values of the intrinsic momentum kt, which is similar to the 
UGDF obtained in [2, 5] at large kt. The inclusion of the hard pp collision within 
LO QCD allows us to describe rather satisfactorily the data at higher values of pt- 
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Figure 3: The inclusive spectrum of charged hadron as a function of pt (GeV/c) in 
the central rapidity region {y = 0) at a/s =7 TeV and wide region of pt (top) and 
the same spectrum at pt < 1.6 GeV/c (bottom) compared with the CMS [14] which 
are very close to the ATLAS data.. 
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Figure 4: The unintegrated gluon distribution xg{x,kt,Qo)/Co as a function of kt 
at x = xo and Qo = l.GeV/c. The sohd hue is the GBW UGDF [2, 5], the dashed 
curve (UGDFl) corresponds to Eq.(ll) and the dotted hne (UGDF2) corresponds 
to Eq.(12) which is very close to the GBW UGDF at h > 1.5 GeV/c. 
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